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Abstract 

The leading-order hadronic contribution to the muon magnetic anomaly = {g M - 2)/2, calculated using a dis- 
persion integral of e + e~ annihilation data and t decay data, is briefly reviewed. This contribution has the largest 
uncertainty to the predicted value of a^, which differs from the experimental value by ~ 3.6 (2.4) standard deviations 
for the e + e~ (t) based analysis. New results since the last workshop and main open issues on the subject are discussed. 
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1. Introduction 

The Standard Model (SM) has been extremely suc- 
cessful. The only missing particle of the SM, the Higgs 
boson, may have been discovered recently at the LHC, 
once verified with more data. All SM predictions have 
been tested often to an extraordinary precision and no 
sign of new physics has been found with few excep- 
tions. One such exception is the well known muon g — 2 
anomaly, a^. The status as of the Tau 2010 workshop 
is about 3.6 standard deviations between the direct mea- 
surement dominated by the E821 experiment at BNL |Q] 
and the corresponding SM predictions J2|- 

The SM prediction a™ is usually decomposed into 
three parts 



fl SM =fl QED +fl weak +a had ; 



(1) 



corresponding to QED, weak and hadronic loop contri- 
butions, respectively. The dominant QED contribution 
includes all photonic and leptonic (e,//,r) loops start- 
ing with the classic a /In Schwinger contribution. It has 
been computed recently through 5 loops and has the fol- 
lowing numerical value |3|: 



af° = (11658471.8951 ±0.0080) x 10 



■10 



(2) 



The weak part includes loop contributions involving 
heavy W ± , Z and Higgs particles. It is suppressed by at 



least a factor ajn ■ m^/M^ 4 x 10 9 . The numerical 
value accounting for the dominant 1- and 2-loop contri- 
butions El |5] 0E] is 



a™ eaK = (15.4±0.1±0.2)xlO 



-10 



(3) 



where the uncertainties stem from quark triangle loops 
and the assumed Higgs mass range between 100 and 
500 GeV, which may be reduced based on the prelim- 
inary Higgs mass determination at the LHC. 

The hadronic part involving quark and gluon loop 
contributions may be further decomposed into leading- 
order (LO), higher-order (HO) and light-by-light (LBL) 
scattering contributions a^ d = a? d ' h0 + aJf d,H0 + 
^had.LBL ^ present, the LO contribution cannot reli- 
ably be calculated from perturbative QCD (pQCD) and 
is determined instead by a dispersion relation |8 1 



.had, LO 



3 (J C 



ds 



(4) 



where R^\s) represents the ratio of the bare cross sec- 
tions of e + e~ annihilation into hadrons to the point-like 
muon-pair cross section and K{s) ~ 1/s is a QED ker- 
nel function |9] and gives a strong weight to low-energy 
part of the integrand. The precision of a^ d ' L0 depends 
thus on that of the e + e~ annihilation data in particular 
that of p(770) — > n + n~ and it has the largest uncertainty 
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to a* and this is why most of the effort from both ex- 
perimental and theoretical sides went into its improved 
determination over the last 20 years or so. 

In the following, we shall briefly describe the new 
development since the last workshop and discuss a few 
open issues on the subject. 



2. New development and open issues 

The preliminary DHMZ 10 results shown at the Tau 
2010 workshop have been published (this and all fol- 
lowing numbers are given in units of 10~ 10 ) ifTOll : 

a had,LO = 692 3 + i.4 + 3.i + 2.4 + 0.2 ± 0.3 (5) 

where the first error is statistical, the second channel- 
specific systematic, the third common systematic, cor- 
related between at least two exclusive channels, and 
fourth and fifth errors stand for the narrow resonance 
and QCD uncertainties, respectively. For this new e + e~ 
based prediction, we included new n + n~ cross section 
data from KLOE, all available multi-hadron data from 
BABAR, a reestimation of missing low-energy contri- 
butions using results on cross sections and process dy- 
namics from BABAR, a reevaluation of all experimen- 
tal contributions using the software package HVPTools 
together with a reanalysis of inter-experiment and inter- 
channel correlations, and a reevaluation of the contin- 
uum contributions from pQCD at four loops. The new 
result is 3.2 below the previous one ifTTI . This shift is 
composed of -0.7 from the inclusion of the new, large 
photon angle data from KLOE, +0.4 from the use of 
preliminary BABAR data in the e + e~ — > n + n^27: mode, 
-2.4 from the new high-multiplicity exclusive channels, 
the reestimate of unknown channels, and the new res- 
onance treatment, -0.5 from mainly the 4-loop term in 
the QCD prediction of the hadronic cross section as well 
as smaller other differences. There was a minor update 
for the FF 2012 workshop [12 | by including the latest 
BABAR 2n + 2n-, 2K2n and 2K2n° channels resulting 
in a£ ad,L0 = 692.4 ± 1.3 ± 3.1 ± 2.3 ± 0.2 ± 0.3. 

The 7i + n~ channel used to be limited in precision, so it 
was proposed in lfT3l to transform the corresponding tau 
spectral function through an isospin rotation to the e + <?~ 
cross section by <x' =I (e + <?~ — > ji~jt~) = Ana 1 / s ■ v(t~ — > 
n~n°v T ) and to provide an independent evaluation after 
accounting for all isospin breaking effects lfl4l . Simi- 
lar transformations can be made for four-pion channels. 
The resulting tau based prediction reads 

fl had, lo [t] = 7Q1 5 + 3.5 + i 9 + 2.4 ± 0.2 ± 0.3(6) 

where the first error is r experimental, the second the 
uncertainty of isospin-breaking corrections iTPfl . the 



third <? + e~ experimental, and the last two the narrow res- 
onance and QCD uncertainties. The 2n and An channels 
account for about 78% of the LO hadronic contribution, 
the rest is taken from the e + e~ channels or pQCD calcu- 
lations. 

Adding to these results the contributions from 
a had,HO _ _ 9 g4 ± 0? com p Ut ed using a similar 

dispersion relation approach, aJj ad LBL = 10.5 + 2.6 ifTBI . 
estimated from theoretical model calculations, as well 
as a^ ED and a™ eak , one gets 



4 M [e + e-] 

_SM 



11 659 180.2 ±4.9,, 



^ M [ T ] = 11 659 189.4 ±5.4,, 



(7) 
(8) 



The e + e (r) based prediction deviates from the direct 
experimental average HI of 



cf^ = 11 659 208.9 ± 5.4 stat ± 3.3 syst 



(9) 



by 28.7 ± 8.0 (19.5 ± 8.3), i.e. 3.6 cr (2.4 cr). 

A compilation of recent a^ M predictions in compar- 
ison with the experimental average of direct measure- 
ments is shown in Fig.[T] In particular the prediction of 
HLMNT 1 1 Ea is similar to that of DHMZ 10 (Eq.((5]l). 
The input e + e~ data sets used are largely identical. They 
differ mainly in the data combination and error treat- 
ment. This is reflected in Table [T] (extracted from Ta- 
ble 4 in ifTBll ). The difference is comparable to or larger 
than one of the quoted errors. In addition, the quoted 
errors are quite different. It is desirable that these dif- 
ferences can be understood and reduced in the future. 



HMNT 07 (eV) 
-287 ±53 

J 07 (eV) 
-294 ±65 

DHMZ 10 (x) 

-195 ±54 

DHMZ 10 (eV) 
-287 ±49 

JS 1 1 (eV+t) 
-292 ± 60 

HLMNT 11 (eV) 
-261 ±49 

BNL-E821 (average) 
0±63 
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Figure 1: Compilation of recent results for (in units of 10 ), 
subtracted by the central value of the experimental average. The 
shaded vertical band indicates the experimental error. 



The difference of 9.1 ± 5.0, i.e. 1.8 cr between the t 
and e + e~ based predictions shown in Eqs.|7]i and ([8}, is 
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Channel HLMNT 1 1 
K + K 22.09 ± 0.46 
n + n- 505.65 ± 3.09 



n n n 



.0 



47.38 ± 0.99 



DHMZ10 
21.63 + 0.73 
507.80 + 2.84 
46.00 ± 1.48 



diff. 
0.46 
-2.15 
1.38 



Table 1 : Comparison for hadronic contributions to in the energy 
range from 0.305 to 1.8 GeV from three K + K~ , n + n~ and 7r + 7!~n° 
channels, extracted from Table 4 in 1151 . 



one of the open issues. Jegerlehner and Szafron claim 
that the difference can be explained by the p° - y mix- 
ing missing in the t data [ 17 1. It remains to be checked 
whether this is the real explanation or there are experi- 
mental issues related to the e + <?~ and t measurements. 
Indeed, the <? + e~ and t difference can be seen from the 
relative shape comparison in the energy range between 
0.3 and 1.4GeV in Fig.[2] 




1.4 

\ls [GeV] 

Figure 2: Relative shape comparison between ALEPH-Belle-CLEO- 
OPAL combined r (dark shaded) and e + e~ spectral function (light 
shaded). 



The other related issue is the different shape between 
BABAR and KLOE n + n~ cross section data (Fig.[3]). 
This difference, leading to an amplified uncertainty in 
the combination following the PDG prescription, pre- 
vents further error reduction. The published KLOE 
measurements were still performed without involving 
the ratio of pion-to-muon pairs as BABAR did. It is 
known that some of the systematic uncertainties cancel 
in the latter ratio measurement. 

Another problematic channel concerns <? + e~ — > 
-2n° (Fig.|4j. There is a large scattering between 
measurements from different experiments, in particular 
between ND and other experiments. In addition when 
comparing the e + e~ average with the r average, there is 
a significant difference in normalization. This discrep- 
ancy deserves further studies and clarification. 
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Figure 3: Comparison between individual e + e~ — » 
tion measurements from BABAR (top) and KLOE 
HVPTools average. 
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3. Running a(s) at M z 

The running electromagnetic fine structure constant, 
a(s) = a(0)/(l - Ai eo (s) - A had (s)), at s = Mf, is an im- 
portant ingredient of the SM fit to electroweak precision 
data at the Z pole. Similar to a^, the error on a(M^) is 
dominated by hadronic vacuum polarization. 

The sum of all the hadronic contributions gives for 
the e + e~ based prediction [ 10 1: 



Aa had (Ml) = (275.0 ± 1.0) x 10" 



(10) 



which is, contrary to the evaluation of aj] ad ' LO , not dom- 
inated by the uncertainty in the low energy data, but by 
contributions from all energy regions, where both ex- 
perimental and theoretical errors are of similar magni- 
tude. This is to be compared with a recent update by 
HLMNT El: Aa had (M|) = (275.5 ± 1.4) x 10~ 4 . 
The reduced electromagnetic coupling strength at Mz 



obtained in Eq.( 10 1 leads to an increase by 7 GeV in the 



central value of the Higgs boson mass obtained by the 
standard Gfitter fit |[T8l to electroweak precision data, 
compared to the previous determination. 
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Figure 4: Cross section of e + e~ — > 2n + 2n° versus center-of-mass 
energy. The shaded green (blue) band gives the HVPTools average 
for e + e~ (r) data. 



4. Summary and perspectives 



The deviation of about 3.6 <x between the direct mea- 
surement and the SM predictions on is significant but 
not sufficient for claiming new physics. The devia- 
tion and the large H — » yy rate observed at the LHC can 
however be explained by a light stau contribution [ 19 1. 

We have mentioned a few open issues in the current 
e + e~ data and the comparison between the e + e~ and r 
data, in particular in the n + n~ and n + n~2ji° channels. 
The tt + ji~ discrepancy between BABAR and KLOE in 
some of the energy ranges prevents us from achieving a 
better precision in the data combination. In order to sig- 
nificantly improve the uncertainty of the leading-order 
hadronic contribution, these issues need to be resolved 
by either more precise new measurements or better the- 
oretical understandings. Lattice calculations are making 
significant progress, but are not yet competitive with the 
dispersion approach with data [20]. The uncertainty of 
the light-by-light scattering contribution is the next item 
to improve. 

The uncertainty of the direct measurement (domi- 
nated by the statistical precision) is now larger than the 
total uncertainty of the SM predictions. Two new g-2 
experiments from Fermilab and JPARC are being built 
and an error reduction by a factor of 4 is expected from 
these experiments in a few years from now. It will be 
a challenge for new SM predictions to match this new 
level of accuracy. 

I am grateful to the fruitful collaboration with my col- 
leagues and friends Michel Davier, Andreas Hoecker and Bog- 
dan Malaescu. 
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